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Abstract 

iV-fold Backlund transformation for the Davey-Stewartson equation is constructed by 
using the analytic structure of the Lax eigenfunction in the complex eigenvalue plane. 
Explicit formulae can be obtained for a specified value of N . Lastly it is shown how 
generalized soliton solutions are generated from the trivial ones. 

Introduction: Inverse scattering transform holds a central place in the analysis of non- 
linear integrable system in either (1+1)- or (2+l)-dimensions [1]. On the other hand it 
has been found that explicit soliton solutions can also be obtained by the use of Backlund 
transformations in a much easier way [2] . These Backlund transformations are also useful 
in proving the superposition formulae for these solutions. There have been many attempts 
to construct explicit A^-soliton solutions for nonlinear integrable system either by Backlund 
transformations or the inverse scattering method. A separate and elegant approach was 
developed by Zakharov et al [3] which relied on the pole structure of the Lax eigenfunction 
and use of projection operators. In this letter we have used an approach similar to that of 
Zakharov et al but have generated a formulae for the A-fold BT of the nonlinear field vari- 
ables occurring in the (2+l)-dimensional Davey-Stewartson equation [4]. Our approach 
is very similar to that of gauge transformation repeatedly applied to any particular seed 
solution. Lastly we demonstrate how non-trivial solutions are generated by starting with 
known trivial ones. 

Formulation: The Davey-Stewartson equation under consideration can be written as 
in + r xx - r yy + r(A 2 - At) = 0, 

iQt + %y ~ Qxx + q(Ai - A 2 ) = 0, ^ 

1 1 

A\ x = - - (q y r + r y q) , A 2y = - - (q x r + r x q) . 
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Equation (|l|) is known to be a result of the consistency of the operators T\ and T 2 written 
as [Ti,T2]^ = 0, where 



2 [ y 1 x yJ \r y / J 2 

To proceed further we set 

M. = fiy j M. = hx 

and 

* = $ exp [i (a + A" 2 ) x-i{(3- A" 2 ) y) . 
Whence the Lax pair becomes, 

M$ = U$, § t = V® (2) 

with 

M — (® x M U-(~ l{a + X ~ 2) ~ q/2 ) (3) 
V0 <W' { -r/2 z(/3-A- 2 ) J (3) 

and 

y _ /^ A + Q( d x> d v) + i hy iQx \ ^ 

\ ik + Q(d x ,d y ) + if 2y J 

With 

A = K 2 - (a + A -2 ) 2 — {p — A -2 ) 2 , 

Q(^, £> y ) = i (d 2 + d 2 y ) - 2 { (a + A" 2 ) ^ - (/? - A" 2 ) 0„} , 

we can construct particular Jost solutions, corresponding to q = qo = const, r = r$ = 
const, A\ = A\q = f® and A 2 = A 2 q = f 2x with /° = f 2x = const. This particular 
eigenvector <3?q turns out to be 



$0 

with 



exp(0ia; + xiy + 6*) exp(0 2 z + X2?/ + &*) 
m exp(6*ix + xiy + n exp(# 2 x + X2y + 6*) 



01 = — i [a + A 2 ) + amo, 6*2 = — % (a + A 2 ) + ano, 
Xl = fe/m + * {P - A" 2 ) , X2 = b/n + i((3 - A -2 ) , 

£i> £2 are arbitrary complex constants, mo, no are arbitrary constants and mo 7^ no- 
Note that det $ 7^ so that ^q 1 exists. 
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Now suppose that <3? n _i denotes the Lax eigenfunction, corresponding to the (n — 
1) soliton solution, and B n (x,y,t) be the transformation which yields the $ n (solution 
corresponding to the n soliton case when applied to <I>„_i), that is 

$ n (x,y,t, A) = B n (x, y, t, A)$ n _i (x, y, t, A). 

Using the above Lax equations @ and (||), we can at once deduce the equations satisfied 
by B„ 

MB n = U n B n — B n U n -i, d t B n = V n B n — B n V n -i, (5) 

where U n , V n denote the Lax matrices. 

Corresponding to the n-soliton solution: Note that U, V are even functions of A, so 
that we can assume that 

B n {-\) = B n {\). 

We now assume B n to have simple pole structure in the complex A-plane, so that 

2A 

B n [x, y, t, A) = Q n + ^ 2 _"^ 2 P n . (6) 
We also assume that 

B- 1 (x,y,t,X)=Q' n + ^^-P^ (7) 

where P n , Q n , P n , Q' n are matrix functions of (x, y, t). The condition B ri B~ l = B~ B n = / 
leads to 

B n (X' n )P^ = 0, P n B~ l (X n ) = 0, 
B- 1 (X n )P n = 0, P' n B n (X' n ) = 0. 

Calculation of the matrices Q n , P n , Q' n , P' n : Let us now go back to equation (|5|) and 
use the expression (|6|) and ([?]). Rewriting U n as 

U n = -iX- 2 I + U' n 

and using the form of B n given in @ we get: 

MQ n + ^^(MP n ) = (-iX~ 2 I + U' n ) (o n + j^P^) 

~ (g» + A^| Pn ) ( " iA " 2/ + U ' n -^ 
which yields equations satisfied by P n and Q n : 

M(P n $„_i(A n )) = U n {X n ){P n $>n-l{Xn)), 

M(Q n $ n _i(A n )) = C/„(A n )(Q n $ n _i(A n )). 

Using the same form of V n as give in equation (^) in the time part, we get the following 
equations: 

9x ( Qn T Pn I — 9y I Q n — P n J , 

d t (P n $„_i(A n )) = U„(A n )(PA-i(A„)), 

Qnt = 0?2Qn + 2i(Qnxd x + Qnydy) + D n Q n — Q n D n ^\, 



352 



S.K. Paul and A. Roy Chowdhury 



where operator OP 2 = i (d 2 + d 2 ) — 2(ad x — (3d y ), a, (3 = const., and D r , 
As per the ansatz of Zakharov we search for P n in the form, 



ifly *9n 

if: 



li- 



ny 



2x 



7nl 
7n2 



It is interesting to observe that 

{S nl ,6 n2 ) = (a„i,a a2 )^l 1 (n), 

where a n ±, a n2 are practically two constants. 
Similarly for P' n and Q' n , we set 



P' 



7„i 

v 7n2 

Whence we get 



(^1X2) and Q' n 



'if 1 + aanSn if 2 + 07722^2 
if 1 + aanhx if 2 + aai 2 5 i2 



where 



F{ 1 = fl\t) exp{imi,(x - y)}, if 2 = /f (i) exp{imi, (x - y)}, 

F? 2 = ff\t)^ V {im' 2l {x-y)}, 



and 



ifW/ (i)exp{zm 2 /Or-y)}, 
2A ; 



07 



2 
AT 



«^if * + « 2 1 

on = — - 



7n 



<7i(<*ji7li +^2^2)' 



^(7u^i + 7/ 2 fe)' 



S12 = 



7/2 



o-K^'i7n + ^2^2) ' 



22 



^(7^/1+7^2^2)' 



12 



P'l + £ Wll 01 + £ Wl2 



2X[ 



HereeJ = — along with 
11 ~ eJ(Mi+W B )' 



' 2 ~ 



7*'i 



7/1 + 7/2 of 



^(7/1^+7/2^) 
Finally the matrix is given as 



72 



7*iA / + 7/2A" 
e J(7/i^i+7/2^ 



12' 



Qi 



/ 2 2 \ 

/z U (*) exp{«m u (x - y)} + —7/1^1 // 12 (0 exp{zm' iz (x - y)} + —7/1^/2 

2 2 

^// 21 (*) exp{zm 2 /(x - y)} + -^7/2^/2 /j 22 (t) exp{im 2 ,(x - y)} + —7^/2 ^ 
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It is also very convenient to rewrite the matrix elements of Q and P in terms of Lax 
eigenfunctions We collect these results below without giving the detailed derivation, 

W| a^ i(A0 _ 6 ^2i i(A0 > Vi a ^22 i(A0 _^2i iW . 



^2 ~~ „ _^22 f\\ u^21 f\\ ' **l 



where 

Af6;/; n ft)exp{imi;(x- j/)} 
H = Af ' 

^ = 4 { AjVj U (t) exp{imi, (a; - y)} + (Af - A?) 6,/, 12 (t) exp^m'^ - y)}} , 

Mi = ± {- (Af - Af ) a,/, n (*) exp{im u (x - y)} - AjV, 12 (t) exp^m'^x - y)}} , 

, Xf ai fl 2 {t) expjim'^x - y)} ^?hff l {t) expjim^x - y)} 
1 ~ Af ' ~ Af ' 

L 'i = jn i^ifiHt) expiimxix - y)} + (Af - Af ) hff 2 (t) exp{zm' a (x - y)}} , 

N l = Ti2 {- i>H ~ Af ) a t ff\t) exp{im 2l (x - y)} - A 2 6,/ 22 (t) exp{zm 2i (x - y)}} , 

AT , _ Afa;/f 2 (t) exp{im 2 i(x - y)} 
l ~ Af • 

The elements of the Pi matrix are: 



pii _ ~ 1 t— iv -i p 

l - ~~ I ^.22 7~\ \ T^2TTVT' i 



^ ^fii(Ai) Dl2 *fii(A|) 



p21 _ ^ ( i— j v p: 

i — ~~ „ ^22 / \ \ u ^21 / \ \ ' i 



a; aj^A,) - b^fi^A,) ' ' a, (^(A,) _ ^(A,) > 
F( $ 22 !(A/) d22 _F/ $2i i( A0 



ai ai^Xi) - bi^XiY 1 ~ a t ai^Xi) - b^XiY 

with 

F l = aifi l {t) ex.p{imu(x - y)} + bifl 2 (t) exp{zm' H (x - y)}, 
F/ = a ; /j 21 (i) exp{im 2i (x - y)} + 6«/f 2 (t) exp{im 2i (x - y)}, 
2Af 



0"; 



^ (Af - A?) 



Construction of the nonlinear fields: Once the form of the matrices Pi and Qi are 
determined we can construct the matrix Bi, so that the Lax eigenfunction ^(A) for the 
next stage can be determined from that of the previous one via, 

$,(A) = B,(A)*,_i(A). 
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These expressions are very complicated, so we just quote one of them to display their 
structure. For example, 



$} 1 (\) = N l /D l , 
A=a,$£ 1 (Ai)-&i$£i(A J ), 

Nt = R^i^i-iW + {R'i - fi(W}*™iM*EiW 

+{M i + /KA)Fj^ 1 (A i )^i 1 ( A ) + M/^ i(A0$ 2i i(A)) 

with similar expression for other elements <1? 12 , <I> 21 and <1> 22 . 
Now, for the determination of nonlinear fields, consider 

U' n Q n = M n Q n + Q n U' n _ 1 , 

where 

In the expression for U' n and U' n _i we take a = (3 = 0, which at once yields 

^Qnx _|_ Qn 

Ql 2 Ql 

This is nothing but a simple recursion relation. Similarly, 



(8) 



7^2 + Til^n-l- (9) 



ggj| Q 22 
Explicitly are can write, 



r » = -# + ® r - 1 - (10) 



2Q{ 2 Q} 1 2Q 2 ^ Qf 

n = 1 > 91 = " of + of «" ri = "oF + oF ro ' 



' ^ 2 r)22 r)22/n22 r)22/Q22 yo ' 

V2 ^2 ^1 ^2 Vl 

2Q| _ 2Q 22 Q 2 ^ + Q22 Q 22 ^ 



So far we have considered f^ff), // 2 (i), f 21 {t), f 22 (t) to be functions of time or constants; 
mu, m 'iii m 2h m 2i to be arbitrary constants; a;, 6/ to be arbitrary constants for all I. Now 
assume that f} 2 (t) = f^(t) = and b t = for all / values. So that Ri = M[ = L t = L\ = 
Ni = F[ = for all I. In this case the form of Bi(X) turns out to be: 



Bi(X) 



J (1 - fef ) fl\t)j m "&-v) (1 + (1 - |r) fl^ty^-^e^i 

f 22 (t)e im2 '( x - y 1 



^i(Xi) 
where <9;_i = 22 . 
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To write the formulae for the n-soliton solution in a compact form, we note that 

{n-R+l I 
i J2 m hn- j+ i{x - y) \ , 

where A%_ R+1 , F£_ R+1 (t) stands for 

n-R+l ^2 n-R+l 

Ai-R+1 = Y\ \k +1 ' K-R+lit) = II /n-i+l(*)- 
1=1 A n-l+l l =1 

Similar expressions can be written for Q„, and its products. Using these, we at once 
obtain: 



KFn{t) exp < i E mi, n - j+ i(x - y) 
-2m n Qf + [_J=i ; 



Gy(t)exp{im' 2in (x - y)} 

G%(t) exp ii }2 m 2 ,n-j+i{x - y) } 

f K [j= ) J 

„_l A \ F K$) exp O £ mi(„_ i+ i)(s - y) > 

El j = l J 19 
i-m n - K Q n _ K , 

K=1 Gr+S) eX P j i E TO 2(n-j+l)(^ - y) > 

where we have set 

if 

m = a(m - n ), G n K (t) = JJ/^+iO), ra^ = -im ln + m! , 

i=i 

m' = &(l/mo - l/n ), ^ n = II/n-mW, T ™ = ^ ( X " ^1) ' 

Qn 2 = T nfl X {t) exp(zmi n x) exp(ra 0:r + m^y + <ft). 
Here <5 is a complex constant. 

Discussions: In the above analysis we have demonstrated how the pole type ansatz of 
Zakharov et al [3] can be used to generate a compact formulae for the iV-fold Backlund 
transformation in the case of the Davey-Stewartson equation. The study yields two main 
results exhibited in equations @ and (|9|) . While the equation gives a recursive procedure 
for the determination of the Lax eigenfunction (<]?; corresponds to the Z-soliton state) 
equation (||) and (|n]) gives the corresponding recursion relation for the nonlinear fields. 
We have actually checked that for n = 1 one obtains the one soliton solution well known 
in the literature. 
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